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Until recently the study of the gravitational field of dark matter was primarily concerned with 
its local effects on the motion of stars on galaxies and galaxy clusters. On the other hand, the 
WMAP experiment has shown that the gravitational field produced by dark matter amplify the 
higher acoustic modes of the CMBR power spectrum, more intensely than the gravitational field of 
baryons. Such wide range of experimental evidences from cosmology to local gravity suggests the 
necessity of a comprehensive analysis of the dark matter gravitational field per se, regardless of any 
other attributes that dark matter may eventually possess. 

In the present note we introduce and apply Nash's theory of perturbative geometry to the study 
of the dark matter gravitational field alone, in a higher-dimensional framework. It is shown that 
the dark matter gravitational perturbations in the early universe can be explained by the extrinsic 
curvature of the standard cosmology. Together with the estimated presence of massive neutrinos, 
such geometric perturbation is compatible not only with the observed power spectrum in the WMAP 
experiment, but also with the most recent data on the accelerated expansion of the universe. 

It is possible that the same structure formation exists locally, such as in the cases of young galaxies 
or in cluster collisions. In most other cases it seems to have ceased, when the extrinsic curvature 
becomes negligible, leading to Einstein's equations in four-dimensions. The slow motion of stars in 
galaxies and the motion of plasma substructures in nearly colliding clusters, are calculated with the 
geodesic equation for a slowly moving object in a gravitational field of arbitrary strength. 

I. DARK MATTER GRAVITY 

The dark matter concept originated from the observed discrepancy between the measured rotation velocity curves 
for stars in a spiral galaxy and clusters and the theoretical prediction from the Newtonian gravitational theory. This 
was first noted by F. Zwicky, when looking at the Coma cluster in 1933 [1]. The measurement of the velocities 
are based on the Tully-Fisher relation between the mass of the galaxy and the width of the 21-cm line of hydrogen 
emissions 0. As fij jure 1 shows, the observed curve becomes almost horizontal (flat), unlike that produced by the 
Newtonian theory [j] . Similar patterns occurs in most spiral galaxies and galaxy clusters [1| . Except for the fact that 
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FIG. 1: Observed Rotation velocity curve of the NGC3198 (error bars) and the prediction from Newtonian theory. 



dark matter interacts mainly through its gravitational field, not much is know about its other physical properties. This 
has motivated several attempts to dismiss dark matter altogether, regarding its gravitational effects as evidences for 
an alternative gravitational theory with respect to Newtonian theory {&, 6], or for variants of general relativity 
General relativity itself has been traditionally excluded from this analysis essentially because the slow motion of the 
observed objects usually leads to the Newtonian limit of the theory. This is reinforced by the fact that near the 
galaxies nuclei where the gravitational field is strong, like that of a black hole, the observed velocities closely agree 
with the prediction from Newton's theory (as in Figure 1.), leading to the conclusion that Newton's gravity should 
apply everywhere else, but a sufficient amount of dark matter must be added to increase the Newtonian gravitational 
pull. A recent review and check list for dark matter candidates can be found in [Til. [T^. 

A qualitative distinction between dark matter gravity and baryons gravity was evidenced in the WMAP experiment 
[l3| . Referring to Figure 2, we quote : ..."Cold dark matter serves as a significant forcing term that amplifies the 
higher acoustic oscillations. Alternative gravity models (e.g., MOND), and all baryons-only models, lack this forcing 
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term so they predict a much lower third peak than is observed by WMAP and small scale CMB experiments"... 
These results were confirmed and sharpened by the observations of the luminous red galaxies in the Sloan Digital 
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FIG. 2: The observed power spectrum of the CBMR of the WMAP experience shows how the dark matter gravitational field can be 
different from the baryonic one. 



Sk^ Survey (SDSS), supporting the AC DAI model [TJ]. This was improved in the recent WMAP fifth year report 

As it seems evident, the study of dark matter gravity and its implications to the formation of structure in the 
early universe must naturally start with gravitational perturbation theory [l6l . [itI [18| . The traditional gravitational 
perturbation mechanisms in relativistic cosmology are plagued by coordinate gauges, mostly inherited from the group 
of diffeomorphisms of general relativity. Fortunately there are some very successful criteria to filter out the latter 
perturbations [l^, [13, |2l[, but they still depend on a choice of a perturbative model. A less known, but far more 
general approach to gravitational perturbation can be derived from a theorem due to John Nash, showing that any 
Riemannian geometry can be generated by a continuous sequence of local infinitesimal increments of a given geometry 

[Hill. 

Nash's theorem solves an old dilemma of Riemannian geometry, namely that the Riemann tensor is not sufficient 
to make a precise statement about the local shape of a geometrical object or a manifold. The simplest example 
is given by a 2-dimensional Riemannian manifold, where the Riemann tensor has only one component -R1212 which 
coincides with the Gaussian curvature. Thus, a flat Riemannian 2-manifold defined by i?i2i2 ~ may be interpreted 
as a plane, a cylinder or a even a helicoid, in the sense of Euclidean geometry. Riemann regarded his concept of 
curvature as deflning an equivalent class of manifolds instead of an specific one [1^. While such equivalence of 
forms is mathematically interesting, it is less than adequate to derive physical conclusions from today's sophisticated 
astronomical observations. 

The solution to Riemann's ambiguity problem was originally proposed in 1873 by L. Schlaefli [25| . conjecturing 
that if a Riemannian manifold could be embedded in another one, then a decision on its real shape could be made 
by comparing the Riemann curvatures of the embedded surface with the one of the embedding space. The formal 
solution of this problem took a long time to be developed and it came only after the derivations of the conditions that 
guarantee the embedding of any Riemannian geometry into another, the well known Gauss-Codazzi-Ricci equations 
of geometry. The most general solution of the Schlaefli's conjecture appeared only in 1956 with Nash's theorem. 

The Gauss-Codazzi-Ricci equations are non-linear and difficult to solve in the general case. Some simplifications 
were obtained by assuming that the metric is analytic in the sense that it is a convergence of a positive power series 
[2^. [27! . Nash's theorem innovated the embedding problem by introducing the notion of differentiable, perturbative 
geometry: Using a continuous sequence of small perturbations of a simpler embedded geometry along the extra 
dimensions, he showed how to construct any other Riemannian manifold [49| [sol. [sH . Nash's approach to geometry 
not only solves the ambiguity problem of the Riemannian curvature, but also gives a prescription on how to construct 
geometrical structures by gradual deformations of simpler ones. 

The purpose of this paper is to show that Nash's geometric perturbative process contributes to explain the formation 
of structures in the universe. To see how this works, in the next section we show how the geometry of the standard 
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cosmology regarded as submanifold embedded a 5-dimensional deSitter bulk, can be perturbed a la Nash, leading to 
a modified Friedman's equation. We draw the theoretical CMBR power spectrum resulting from such perturbation, 
comparing it to the observed spectrum in the WMAP experiment. In section 3 we apply the same procedure for local 
dark matter gravity. However, since there are no evidences that the perturbative process is still going on today, except 
perhaps in young galaxies and cluster collisions where experimental data is still scarce, we assume as a first estimate 
that Nash's perturbation is locally negligible in most spiral galaxies. This leads to a simpler set of equations where the 
brane-world gravitational equations reduce to the usual four-dimensional Einstein's equations. In this case we will see 
that the motion of stars and of plasma substructures can be described by the application of " slow geodesic motion" 
equation which holds true in an gravitational field of arbitrary strength and is not limited to the weak gravitational 
fields as in General Relativity. 

II. DARK MATTER COSMOLOGY 

We start by reviewing some basic ideas of Nash's geometric perturbation theorem: Suppose we have an arbitrarily 
given Riemannian manifold V„ with metric 5^^, which is embedded into a given higher dimensional Riemannian 
manifold Vd, the bulk space. Then we may generate another metric geometry by a small perturbation g^i, + Sg^j^^, 
where fso} 

= -2fc^,,52/^ a = N + l...D (1) 

where Sy"^ denotes an infinitesimal variation of the extra dimensions orthogonal to Vn and fc^^a denote the extrinsic 
curvature components of Vn relative to the extra dimension y". Using this perturbation we obtain new extrinsic 
curvature kp,^a , and by repeating the process we obtain a continuous sequence of perturbations like 

g^Li' = (jfii' + Sy"' k^^a + Sy°-Sy'' g'"^k^pahab ■ ■ ■ 

In this way any Riemannian geometry can be generated. 

Nash's original theorem used a flat D-dimensional Euclidean space but this was soon generalized to any Riemannian 
manifold, including those with non-positive signatures j!!]. Although the theorem could also be generalized to include 
perturbations on arbitrary directions in the bulk, it would make its interpretations more difficult, so that we retain 
Nash's choice of independent orthogonal perturbations. It should be noted that the smoothness of the embedding 
is a primary concern of Nash's theorem. In this respect, the natural choice for the bulk is that its metric satisfy 
the Einstein-Hilbert principle. Indeed, that principle represents a statement on the smoothness of the embedding 
space (the variation of the Ricci scalar is the minimum possible). Admitting that the perturbations are smooth 
(differentiable), then the embedded geometry will be also differentiable. 

The Einstein-Hilbert principle leads to the D-dimensional Einstein's equations for the bulk metric Gab in arbitrary 
coordinates 

T^AB - -^TLQab = oi*TXB (2) 

where we have dispensed with bulk cosmological constant and where TJ^ denotes the energy-momentum tensor of 
the known matter and gauge fields. The constant a* determines the D-dimensional energy scale. 

The four-dimensionality of the space-time manifold is an experimentally established fact, associated with the Poincar 
invariance of Maxwell's equations and their dualities, later extended to all gauge fields. Therefore, all matter which 
interacts with these gauge fields must for consistency be also defined in the four-dimensional space-times. On the 
other hand, in spite of all efforts made so far, the gravitational interaction has failed to fit into a similar gauge scheme, 
so that the gravitational field does not necessarily have the same four-dimensional limitations and it can access the 
extra dimensions in accordance with ([T]), regardless the location of its sources. 

We assume that the four-dimensionality of gauge fields and ordinary matter applies to all perturbed space-times, 
so that it corresponds to a confinement condition. In order to recover Einstein's gravity by reversing the embedding, 
the confinement of ordinary matter and gauge fields implies that the tangent components of a^TJ^ in the above 
equations must coincide with SttGT^^ where T^i^ is the energy-momentum tensor of the confined sources [5l| . 

Since dark matter gravity is not necessarily confined, it also propagates in the bulk. Furthermore, based on the lack 
of experimental evidences on the energy-momentum of the dark matter in the bulk, we also assume that the normal 
and the cross components of the dark matter energy-momentum tensor, respectively T^a and Tab vanish, meaning 
that there are no known sources outside the four-dimensional space-times. 
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The standard Friedman-Lemaitre-Robertson- Walker (FLRW) universe can be embedded without restrictions in 
a five-dimensional. Therefore, using the Einstein equation ^ written in the Gaussian frame defined by the four- 
dimensional submanifold, we obtain the equations of the FLRW embedded geometry [30l| 

Rfii^ — -^Rgt^u — Qtj.u = —SttGT^,, (3) 

fc^;p-^M = (4) 

where now T^i, is the energy-momentum tensor of the confined perfect fluid, fc^i^s = k^^^, denotes the components of 
the extrinsic curvature of the embedded space-time, h = g^'^k,y^, = k^^'^k^^ and 

Q^j.v — k^fj^kpi, hk^n 2^"^ ^ (^) 

This tensor is independently conserved, as it can be directly verified that (semicolon denoting covariant derivative 
with respect to g^u) 

Q^",. = (6) 

In coordinates (r, 0, t) the FLRW model can be expressed s [43|: 

= -dt^ -t- a2(t)[dr2 4- fi{r){d0^ + sm^ed(f>^)] (7) 

where /fe(r) = r, sinr, sinhr corresponding to fc = 0, -f 1, —1 (spatially flat, closed, open respectively). We start solving 
(m for the above metric in the deSitter bulk. It is easier to flnd the solution of Codazzi's equations 

of which (HI) is just its trace. The general solution of this equation is 

= ^44 = -T^f-) i,j = 1...3 (8) 



Defining B ~ -r, we may express the components of Q^u as 



u B \ 36*^ (yb^^ B 

<3u= ^ ( 2— - ij gy, Qi4 = --^, Q = --^J^^ i,j = 1..3 (9) 



where H = a/a is the usual Hubble parameter. After replacing in ^ we obtain Friedman's equation modified by the 
extrinsic curvature: 



,a 2 k SvrG 6^ 
a a-^ 3 



To interpret this result we have compared it with the XCDM, phenomenological x-fluid model, with state equation 
Px = oJxPx, which corresponds to the geometric equation on b{t) 

I^hi-Sivx)- (11) 

2 a 

This cannot be readily integrated because Ux is not known. However, in the particular case when ujx = <j;o=constant, 
we obtain a simple solution 

b{t) = 60(^)^(1-3^0) (12) 
ao 

where oq and &o 7^ are integration constants. Replacing this solution in (|10p we obtain an accelerated universe which 
is consistent with the most recent observations, when the values of loq are taken within the range —1 < ujq < —1/3 
[30| . Furthermore, the theoretical power spectrum obtained from the extrinsic curvature perturbation of the FLRW 
model, is not very different from the observed power spectrum in the WMAP/5y, displayed in Fig. 2. 
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FIG. 3: The theoretical power spectrum calculated with the CAMB for —1 < luq < —1/3, Massive Neutrinos=l, masslcss neutrinos =3.04. 



III. LOCAL DARK MATTER GRAVITY 

It is possible that in young galaxies which are still in the process of formation [3^ ; in galaxies with active galaxy 
nuclei [33|; or even in cluster collisions, Nash's perturbations could be applied, where the metric symmetry is taken 
to be local. However, there not sufficient experimental data to support such local perturbations. Therefore, in the 
following we restrict our study on local dark matter gravity to the cases where Nash's perturbations has ceased. From 
^ we see that this limiting condition is given by k^^^. Furthermore, if T^i, is composed by ordinary matter and gauge 
fields, the rsulting equations are the same as the vacuum Einstein's equations 



R^. - 7:R9^. = (13) 



These equations are be understood in the context of the embedded space-times and with the confinement conditions 
for ordinary matter and gauge fields. They do not represent the whole of General Relativity because the princi- 
ple of general covariance does not necessarily apply to the bulk geometry. This follows from the fact that Nash's 
perturbations are restricted to be along the orthogonal directions only. 

The recently observed correlation between the highest-energy cosmic rays sources and active galactic nuclei [s^ 
suggests the existence of strong gravitational fields in those regions [35. 36, 37). Nonetheless, the observed motion 
of stars near black holes observed at the nuclei of some of these galaxies indicate that the velocities are small, of 
the order of a few hundreds of kilometers per second. This requires the description of a slow geodesic motion in the 
presence of strong gravitational fields In what follows we use essentially the first part of the exposition in the Misner, 
Thorn & Wheeler's book 

Consider a slow free falling particle (or star) in a gravitational field described by a solution of the vacuum Einstein's 
equations (like in (|13p in a suitable coordinate system. Initially the particle is located in a far away region, where 
the action of the gravitational field is supposedly weak, and we can write (The notation /i^^ reminds that this is not 
Nash's perturbation) 

+ Shfj_^, where << (14) 

Here the small deviation from Minkowski's metric has nothing to do with the velocity of the particle. However, since 
V << 1, we can use Newtonian coordinates (with x'^ = t, t being the Newtonian time), so that the spatial components 
of the geodesic equations become 



dr'^ dr" rfr^ 1 



dt^ dt dt dt 2 

At this point a scalar field Lp may be defined such that 



dt'^ dx^ 

Comparing (|15p and (|16p we obtain an equation to determine the field (p: 



(16) 



dip 1 dShii 
9^ = -2^'^^ 
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As the particle continues its fall, the gravitational pull continuously builds up by small increments of the metric as in 

Qtiu ~ Vt^i' + Sh^,u + {Sh^^f H 

Actually, there is no way to stop this process without applying an external force. Thus, equation (jl7p can be 
integrated along the geodesic, where Sh^i, continuously increase up to a finite value h^^, leading to the scalar gravi- 
tational potential (often referred to as the nearly Newtonian potential, not to be confused with the post Newtonian 
approximations) : 

^^~2Jo ^(^^^44) = -2(1+544) (18) 

Notice that except at the beginning of the free fall a weak gravitational field was not imposed. Here 544 is obtained 
from an exact solution of Einstein's equations, so that at the end all metric components 17^1/ contribute to (|18p . In 
the following we exemplify this application of to the motion of stars in galaxies and clusters. 
%subsubsectionRotation Velocity curves in Galaxies 

The gravitational field of a simple disk galaxy model can be obtained from the cylindrically symmetric Weyl metric, 
expressed in cylindrical coordinates {r,z,9,t) as [39j : 

where A = A(r, z) and a = a{r, z). As shown in [40l . [4l| . the Weyl cylindrically symmetric solution is diffeomorphic 
to the Schwarzschild solution. Replacing the Weyl metric with these conditions in Q and ([6]) we obtain 

-A.^ + rcT^ - ra^^ = (20) 

— (T_r — ra^rr — fO^zz — (21) 

A,„ + A,,, + (T^, + (J%=0 (22) 
2ra,rO,z — A_2 (23) 

The cylinder solution is diffeomorphic to a Schwarzschild's solution [s^. 

In the following we apply this solution to find a the geodesic motion of slowly free falling star under the gravitational 
field in the galactic plane. Since the slow motion geodesic equation (jl6p is not invariant under diffeomorphisms, we 
may consider three separate stages separately, in accordance with the symmetry of the gravitational field which is 
effective at the current position of the star: 

(A) The star is far away from the galaxy: 
In this case, the gravitational field of the galaxy is weak, like that of a point source. Therefore, the predominant 
gravitational field is given by the exterior Schwarzschild solution, as seen from a large distance from the galaxy 
nucleus. In spherical coordinates, we can write the metric component (744 = —(1 — 2Af/r), and (|18p gives 

v\Newton = -M/r, r > > Tq 

which is equivalent to the Newtonian gravitational potential produced by a distant mass M, determined by the New- 
tonian (weak field) limit. In such situation (jlSp agrees with all estimates resulting from the Newtonian gravitational 
theory for a spherically symmetric dark matter halo, producing the same rotation velocity curves. 
(B) The star is close to the galaxy's disk: 

When the star is close to the galaxy, in the disk plane, the predominant dark matter gravitational field can be sim- 
ulated by the vacuum cylindrically symmetric Weyl metric, satisfying the vacuum Einstein's equations ([^0]) to (|23p . 
However, to configure a disk galaxy we need to apply the condition that Weyl's cylinder has a thickness that is much 
smaller than its radius: h « rg. With this condition we can no longer apply the diffeomorphism invariance of general 
relativity, so that the solution must be written in cylindrical coordinates. 

The disk-symmetry condition can be written as z € [— /io/2, /io/2], for r e [0, tq], /iq << fo, so that the 
functions a{r,z) and \{r,z) may be expanded around z = as 

cr(r, z) = a{r, 0) + za{r) + ■ ■ ■ 
A(r, z) = A(r, 0) + zb{r) H 

where we have denoted a{r) — ^'^j^'^^ and b{r) — ^^j^^''^ . Neglecting the higher order terms, it follows 

that equations (PT|) and become a simple system of equations on cr, with general solution (T(r, z) = Inr + ci{z) 
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where 02(2) is an r-integration constant and where we have denoted K = b{r)/a{r). Derivation of a with respect to 
z gives C2{z) = a{r)z + cq, but since C2 does not depend on r, it follows that a(r) must be a constant Qq. By similar 
arguments we find that b{r) = bo =constant, so that K — Kq ~ &o/ao is also a constant. Replacing these results in 
(PO)) and (im, we again obtain another simple solvable system of equations in A, so that the solution of the vacuum 
Einstein's equations for the Weyl disk is 

a{r,z) = ^\nr + aoz + Co (24) 

X{r,z) = ^\nr~aoY + boz + do (25) 

where Cq and do are again integration constants. From ([24|) we obtain 1744 = —e^'^ = —g^^ '"''e^^o^e^'^''. Therefore, 
for a star near the galaxy in the galaxy plane 2; = 0, in the region between the nucleus radius Vc and the disk radius 
ro (Ull) is 

<^Jd^sfe--i(l+g44)J.=o = -^(l-e2"V^''), rc<r<ro (26) 

Here we cannot make use of the Newtonian limit to determine the integration constant e^'^" in ([26|) . because we do 
not have the same symmetry and the same boundary conditions appropriate for the Newtonian gravitational field. 
Instead, we may compare (|26)) with the local Newtonian potential produced by a disk of visible mass M, suggesting 
that the above integration constant can be written as proportional to the to the total baryonic mass Mh of the galaxy. 
Thus, in units G=c=l we set e^'^°Ko/2 — f3oMb, where f3o is a mass scale factor, proportional to the estimated visible 
mass of the galaxy [43j] . 

The rotation velocity of a test particle under is v = ujor, wq = constant, obtained by comparing the radial 
force ^r, with —^f, so that the velocity is given by w = \fV^\- the particular case of the disk, using (pS]) we 
obtain 

v{r) = ^J\|3oMbr^^>l < r < ro (27) 

Here Kq represents the ratio between the coefficients a(r) and b{r) of the expansion of the metric functions a and A 
respectively. In the considered linear expansion in and (^5]) . Kq can be adjusted experimentally. In practice it 
can be scaled by a constant, so as to belong to interval between and 1, corresponding to the galaxy nucleus radius 
where the spherical symmetry applies, and to the distance where the disk symmetry applies. The extreme values in 
that interval are excluded: The value Kq = is excluded because for that value the potential (pS)) does not exert any 
force. The value Ko = 1 is also excluded because it gives velocity proportional to ^Jr which is not experimentally 
verified. 

Using the minimum square root curve fitting, figures 2, 3 and 4 show the velocities near the disk, calculated with 
([27|) for some known galaxies (red dots) for r larger than the estimated core radius where the disk symmetry applies. 
The values of the two free parameters Ko and (3o in each case were determined by the minimum squares numerical 
iteration method, extracted from known experimental data [45j . For comparison purposes, we included the black 
error bars showing the measured velocities and the green triangles corresponding to the Newtonian predictions. 

(C) The star is near the galaxy's nucleus: 
When the star reach the nucleus of the galaxy, the observed velocity is still small, but the gravitational field is likely 
to be strong. A simple configuration of the gravitational field acting on the star can be obtained if we neglect the 
rotation of the galaxy's nucleus, so that we obtain a spherically symmetric gravitational field, described by the by the 
Schwarzschild solution of (fT^l) with 1744 = — (1 — 2M/r). Then we obtain from ([T5)) 



J Nucleus 



^^-M/r, r<rc (28) 



and the star rotation velocity is given by v(r) — ^ M jr^ were now M is associated with the Schwarzschild mass 
of the Nucleus. Thus, the velocity looks exactly like the one in the Newtonian theory, with the exception that the 
gravitational field is not necessarily weak. This is a consequence of the assumed metric symmetry at the galaxy's 
nucleus. 

Two recent observations of colliding cluster dark matter halos have provided new insights on the dark matter gravity 
issue: One of them is the bullet cluster 1E0657-558, showing the motion of a sonic boom effect in the intercluster 
plasma with velocity ~ 4700km/s, visible through x-ray astronomy (46l |. Under the assumption that the gravitational 
field of the two clusters is Newtonian and the existence of gravitational field of dark matter halos in each cluster, 
the (linear) superposition of the two dark matter gravitational fields produce a center of mass of the system which 
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FIG. 4: Via Lactea with Mb = 1.99 x 10"*' Kg, R = 8, 5Kpc, Ko = 0.0682 and f3o = 0.0901. 
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FIG. 5: NGC3198, with Mb = 1.19 x 10^"Kg, R = 16, 6Kpc, Kq = 0.0162 and (3q = 0.8224 
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FIG. 6; NGC3949, with Mb = 4.97 x IQ-'^Kg, R = 6Kpc, Ko = 0.3766 and /3o = 1.1657 



coincides with the observed position of the plasma bullet. This has been claimed to be an observational evidence for 
the existence of dark matter. 

On the other hand, using (|18p instead of the Newtonian gravity, and admitting that the gravitational field of each 
cluster is spherically symmetric, we have the gravitational field of the spherically symmetric dark matter halos would 
be given by two separate Schwarzschild solutions. Then the action of each of these gravitational field on the plasma 
can be taken as in ([28|) . which resembles the Newtonian potential, producing the same center of mass as described 
in Newtonian theory. However, such superposition of two Schwarzschild's solution is only a crude approximation to 
compare with the dark matter halos. To be more precise, either we consider a two body problem solution of (jl3p . 
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or else we consider the cluster collision as a Nash's perturbation process. In the latter case we may start with an 
embedded spherically symmetric cluster which is perturbed by the extrinsic curvature generated by the second cluster 
according to This is easier said than done, because the differentiable embedding of a Schwarzschild solution 
requires six dimensions. Nonetheless, in principle this solution can be calculated and tested. 

The second observation is that of the Abel 520 cluster MS0451+02, showing again two colliding "dark matter halos" , 
where at least one of them do not seem to be anchored to a baryonic structure. The most immediate explanation is 
that this "pure dark matter halos" may be an evidence of a non- linear effect of the dark matter gravitational field, 
which does not agree with the Newtonian gravitational field assumption. So, in a sense this observation backs up the 
hypothesis of a non-linear alternative gravitational theory at the galaxy/cluster scale of observations. 

In particular it is possible to explain the existence of a gravitational effect which is not anchored to a baryonic 
source as a solution of ([T^ . as a for example a vacuum Schwarzschild solution or a Weyl disk solution, which may act 
as a perturbation to another cluster, again applying ([T]) to find the final gravitational field of the system. The result 
may be compared with the motion of the x-ray observed plasma structure. However, we are still pending on further 
details on the Abel 520 collision. 



IV. SUMMARY 

Gravitational perturbation theory has become an essential tool for the explanation of the formation of large struc- 
tures in the universe. Here we have presented an application of Nash's theorem on perturbations of geometries to 
the formation of space-time structures and to explain the local effects usually attributed to dark matter. In a brief 
justification we have argued that the theorem improves Riemann's geometry in the sense that it replaces the somewhat 
absolute notion of Riemann curvature, by a relative notion of curvature with respect to the Riemann tensor of the 
bulk defined by the Einstein-Hilbert principle. 

The relevant detail in Nash's theorem, is that it provides a mathematically sound and coordinate gauge free way 
to construct any Riemannian geometry, and in particular any space-time structures, by a continuous sequence of 
infinitesimal perturbations along the extra dimensions of the bulk space, generated by the extrinsic curvature. 

The four-dimensionality of space-time is regarded here as an experimental fact, related to the symmetry properties 
of Maxwell's equations or, more generally of the gauge theories of the standard model. Thus, any matter that interacts 
with gauge fields, including the observers, must remain confined to four-dimensions. However, in accordance with 
Nash's theorem, gravity as represented by a metric cannot be confined because it is perturbable along the extra 
dimensions. 

In a first cosmological application of that theorem we have compared our results with the present observational data. 
Starting with the FLRW cosmological model embedded in the five-dimensional deSitter bulk, and applying equations 
^ and We have found that Friedman's equation is necessarily modified by the presence of the extrinsic curvature. 
We have shown that this modification is consistent with the observed acceleration of the universe, including with the 
observed power spectrum of the CMBR. 

In principle the local effects of dark matter gravitation should have the same explanation, although they have 
different observational aspects. The simplest case is that of the rotation curves in galaxies and galaxy clusters, which 
originated the dark matter issue. In this case, there is no evidence that the gravitational perturbation process is still 
active, except perhaps in young galaxies. Therefore, we have considered the case where Nash's perturbations vanish, 
obtaining the vacuum Einstein's equations, applied to the metric with an specific symmetry. Instead of the Newtonian 
potential we have applied the geodesic equations for slow motion. Using the Weyl metric to simulate a disk galaxy, 
we have compared the result with some known cases. 

More recently the observations of merging clusters have provided additional information on the dark matter issue. 
The slow motion of the plasma substructures can be handled by the same equations but it appears to us that the 
correct formulation of the problem should be made with Nash's perturbative analysis. 
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